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Abstract 

A family of the so-called Maxwell times which arises in the context of Higher-Order Generalized 
Hydrodynamics (also called Mesoscopic Hydro-Thermodynamics) is evidenced. This is done in 
the framework of a HOGH build within a statistical foundation in terms of a Non-Equilibrium 
Statistical Ensemble Formalism. It consists in a description in terms of the densities of particles 
and energy and their fluxes of all orders, with the motion described by a set of coupled nonlinear 
integro-differential equations involving them. These Maxwell Times have a fundamental role in 
determining the type of hydrodynamic motion that the system would display in the given condi¬ 
tion and constraints. The different types of motion are well described by contractions of the full 
description done in terms of a reduced number of fluxes up to a certain order. 
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I. INTRODUCTION 


We present a generalization of Maxwell time [1], which can also be referred-to as gen¬ 
eralized Maxwellian relaxation time [2], a family of which is present in the evolution equa¬ 
tions of Higher-Order Hydrodynamics, also dubbed as Mesoscopic Hydro-Thermodynamics 
(MHT) [3-6]. Such generalized Maxwell times are associated to the dampening of densities 
of particles (molecules in a classical fluid) or quasi-particles (Bloch-band electrons, harmonic 
phonons, magnons, polaritons, etc., in solid state matter, in semiconductors, for example), 
and of their energy density, together with the fluxes of all orders of both. 

The origin of Maxwell times goes back to the fundamental article by J. C. Maxwell in 
1867 on the dynamical theory of gases [1], in what is related to viscoelasticity and which 
is being expressed as the quotient of the dynamical viscosity coefficient with the modulus 
of rigidity. According to Maxwell it may be called the ’’time of relaxation” of the elastic 
force. It has been noticed [2] that given a fluid subjected to some variable external forces, 
which vary periodically in time with frequency co, of the period 1/co is large compared with 
Maxwell time , i.e. cod <C 1, the fluid under consideration will be have as an ordinary viscous 
fluid. If however, the frequency co is sufficiently large (so that cod S> 1, the fluid will be have 
as an amorphous solid. 

We present a detailed study of the generalized Maxwell times that arise in a Meso¬ 
scopic Hydro-Thermodynamics, built upon a Non-Equilibrium Statistical Ensemble Formal¬ 
ism (NESEF for short) [7-10], at the classical and quantum mechanical levels of HOGH 
described in Refs. [5,6,11,12], 

In Section II is considered the classical MHT of a fluid of particles embedded in a thermal 
bath. In Section III the case of a MHT at the quantum level of a system of phonons in 
semiconductors, and in Section IV a MHT at the quantum level of a system of electrons in 
doped semiconductors. 
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II. MAXWELL TIMES IN THE MHT OF A CLASSICAL FLUID 


Construction of MHT in the framework of a Nonequilibrium Statistical Ensemble Formal¬ 
ism is presented in Ref. [5]. There it is provided an approach enabling for the coupling and 
simultaneous treatment of the kinetic and hydrodynamic levels of descriptions. It is based 
on a complete thermo-statistical approach in terms of the densities of matter and energy and 
their fluxes of all orders, as well as on their direct and cross correlations, covering systems 
arbitrarily far-removed from equilibrium. The set of coupled nonlinear integro-differential 
hydrodynamic equations is derived. They are the evolution equations of a Grad-type ap¬ 
proach involving the moments of all orders of the single-particle distribution fi(r, p;f), 
derived from a generalized kinetic equation, built in the framework of a Nonequilibrium 
Statistical Ensemble Formalism [13]. 


The moments of the single-particle distribution function, in momentum space p, are the 
hydrodynamic variables 



n(r,t) = J d 3 .p/i(r,p;t), 

(1) 

which is 

the densities of particles, 



I„(r,f) = J d 3 p u(p)/i(r,p;t), 

(2) 

with 

u(p = — , 

(3) 


m 

where I„ 

is the flux (current) of particles, 



4 2] ( r >*) = J d 3 pu [2] (p)/i(r, p; t ), 

(4) 


where = [u : u] is the inner tensorial product of vectors u, with being the second- 
order flux (or flux of the flux), a rank-2 tensor, which multiplied by the mass is related to 
the pressure tensor, and 

I [ n ] (r,t) = I d 3 pu [f] (p)/i(r, p; t) , (5) 

are the higher-order fluxes of order l ^ 3 (the previous three of Eqs. (1), (2) and (4), are 
those for t = 0,1 and 2 respectively) where is the Crank tensor consisting of the inner 
tensorial product of i vectors u of Eq. (3), that is, 

u [f] (p) = (£— times)(6) 

Lm m m. 
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We do have what can be called the family of hydrodynamical variables describing the 
material motion , i.e., the set 


Mr,*); In(r,t); {/^(r,t)}} , 


(7) 


with £ = 2, 3,..which we call the MHT-family n. 

On the other hand, we do have the family of hydrodynamical variables describing the 
thermal motion , which we call the MHT-family h, consisting of 

.2 


h(r,t) = I d 3 p|^/i(r,p;t) , 

I h {r,t) = [ d 3 p^——f 1 (r,p-,t), 

J 2 m m 

A ?(M) = f ^ 3 p^ uM (p)/i( r , p; *), 

with £ = 2, 3,..that is, in compact form, those in the set 

W r ,*); I h(r,t)-, {/i fl (r,«)}} , 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 


which are, respectively, the density of energy, its first vectorial flux (heat current), and the 
higher-order tensorial fluxes. It can be noticed that in this case of a parabolic type energy 
momentum dispersion relation, E(p) = p 2 /2m, the set of Eq. (11) is encompassed in the 
previous one: in fact 

( 12 ) 

(13) 


MM) = yIY{4 2 l(r, t )}, 


W r »l) = ^rr 2 {4 3, (r,*)} , 

where Tr 2 stands for the contraction of the first two indexes, and, in general 

M) = ?^r 2 {/>f +2] (r,*)} , 


(14) 


for all the other higher-order fluxes of energy, that is, any flux of energy of order £ is contained 
in the flux of matter of order £ + 2. 

Let us consider the equations of evolutions for the basic macrovariables of the family of 
particle motion, which for the general flux of order £ {£ = 0,1, 2,...) is 


J^ ] (r,t) = f d 3 pu M (p)^/i(r,p;f), 


(15) 
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Using in Eq. (15) the evolution equations for /i(r, p; t) of Ref. [13] it follows the general 
set of coupled equations for the density, 7 = 0, the current, £ = 1, and all the other higher- 
order fluxes, £ ^ 2, given by 

7/6 

€ 

a L0 ^(i, S ) ( V/ n _1] ( r ’ <)] + 2^1 V • /Jf +1] (U t) + 
s=l 

+ SP(r,i) ■ (16) 

The last term on the right of Eq. (16) is given by 

= b T0 {a(l, s)[l [2] I^- 2 \r,t)} + b Tl — (a(l, s)[l l2] /^“ 2] (r, f)} 

Tfl 

£ 

+3£a r i—I$(r,t) + a L1 — ^ a(l, s)[V/^" 1] (r, f)] + R [ ^(r,t) , (17) 

777 / 77 7 

5=1 

with all quantities involved given in Ref. [5]: we omit the details that are not necessary for 
the analysis of Maxwell times. 

The Erst contribution on the right of Eq. (17) contain the flux of particles of order 7 — 2, 
whereas the other three contributions are associated to the fluxes of energy of orders 7 — 2, 
7 — 1 and 7, terms that can be considered as associated to thcrmo-striction effects which 
couple these equations with the set of kinetic equations describing the motion. However, it 
can be noticed the already mentioned fact that the fluxes of energy can be given in terms 
of those of particles, as given in Eqs. (12) to (14). 

The second term on the right of Eq. (16), for the evolution of the flux of order 7, is the 
one that contains the Maxwell time 9 n j> associated to such flux, namely 


0 n £ ~ ~ 7 «r0 + 7(7 — l)|& T i| , 

(18) 

with the coefficients a r0 and b Tl given by 


= f dQQ*f r0 (Q ), 

(19) 

with 

, ,„\ n R M$ /2 ir\ip(Q)\ 2 fm ^ 

fAQ) - o U + 1 b 

(20) 
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where i/)(Q) is the Fourier transform of the potential energy w(\rj — R M |), between the j'-th 
particle at position r and the one in the thermal bath at position R /I; , ur is the density 
of particles in the thermal bath, V is the volume, and (3q 1 = A^Tq, with T 0 being the 
temperature of the thermal bath, moreover 


b T 1 


CLtO 

~Y 



( 21 ) 


and m and M are the mass of the particles and of the particles in the thermal bath respec¬ 
tively. According to Eqs. (18) and (21) we can write 


9 


-i 

ni 


^I^tO | 


i + i(«-D 


M V 
m + M ) 


( 22 ) 


and it can be noticed that |a T o| is the reciprocal of the Maxwell time of the Erst flux which is 
the reciprocal of the relaxation time of the linear momentum p(r, t), once p(r, t) = ml n (r, t). 

Maxwell times of Eq. (18) have their origin from the collision integral J?\r, p;f) in the 
kinetic equation for the single-particle distribution present in Eq. (A. 11) in Appendix A in 
Ref. [13]. It arises out of the interaction with the thermal bath in the contributions called 
effective friction force and diffusion in momentum space. 

In a recurrent procedure it follows from Eq. (18) that we can write 




1 + 


M 




^nl i 


(23) 


5 (m + M) 

what tell us that any characteristic time for t ^ 2 is proportional to the one of £ = 1, that 
is, the one for the first flux which, as notice above, multiplied by the mass m is the linear 
momentum density and then all are proportional to the linear momentum relaxation time. 
On the other hand we do have that 


9ne -|-i f 5(1 + x) + £ — 1 

0 ne i + 1 5(1 + x) + i 

for i = 1, 2, 3,... and where x = m/M, then the ordered sequence 


(24) 


9 ni > 0 n „ > 0 


n 2 


'ns 


> o ne > 0, 


ne +1 


> 


(25) 


is verified, and it can be seen that 0 nt —» 0 as i —> oo. Moreover, according to Eq. (23) it 
follows that 


9 


n e 


5(1 +£) „ 

<[5(l + i)+<-l) 


(26) 
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Comparing with the second flux (£ = 2) which is the one related to the pressure tensor, 
it follows that for Brownian particles (x 3> 1) 9 n2 ~ Q n J 2 and for Lorentz particles (x <C 1) 
0 n2 ~ 59 ni /12. A comparison with the third flux leads to the results that 9 n3 ~ 9 ni /3 and 
6 n3 ~ 59 ni /21 for the Brownian and Lorentz particles respectively. For any i we do have 
approximately: 

1) for the Brownian particle ( m/M 3> 1) 


9 


m 



(27) 


2) for the Lorentz particle (m/M <C 1) 

k 

0 rs '~'- o 

nt - (4 + £)l ' 


(28) 


or 9 nt ~ 5 9 ni /l 2 for large £. 

According to Eq. (24) as the order of the flux largely increases its Maxwell characteristic 
time approaches zero, and 9 nt+1 /9 nt ~ 1, with both practically null. In Fig. 1 it is displayed 
the ratio of some characteristic times compared with the one of the linear momentum in 
terms of the ratio m/M. 



FIG. 1: The quotient between Maxwell characteristic times and the one of the first flux as a 
function of x = m/M. 


We emphasized that the hierarchy shown in Eq. (25) is evidencing the important fact of 
the ever decreasing value of the Maxwell times with the increasing order of the fluxes. This 
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is quite relevant for the question of how to introduce a contracted description (i.e., keeping 
a reduced number of fluxes) for characterizing the particular type of hydrodynamic motion 
(diffusion, damped waves, etc.) that the system displays, as shown later on. 



III. MAXWELL TIMES IN THERMAL TRANSPORT 


In this section we concentrate the attention on the evolution of the thermal quantities 

(MM); Mr./); Uj?M)}}* (29) 

with d = 2, 3,..that is, the density of energy and its fluxes of all orders given in Eq. (11). 
Their evolution equations are 

M) = ’ ( 30 ) 

in which we need to introduce the kinetic equation for the single particle /i(r, p;f) which 
is given in Ref. [13]. Performing the lengthy calculations involved we finally arrive at the 
general evolution equations {£ = 0,1, 2,...) 

l 

= -^(r,t)4* +1] (r,t) - J2 a ( 1 i s )[F{r, t )I l h~ 1] (r,t)] + d^lP{r,t) + 

S= 1 

t 

a L o 5] <r(l, s)[ViT 1] (r, t)] + 2(£ + 3 )a L1 V • /J? +1] (r, *) + 

S=1 

Sj?M), (31) 

where 

6hi — (£ + 2)|a r o| + [£ 2 + 5(£ + l)]|6 r i|, (32) 

is the reciprocal of the Maxwell time associated to the f'-th order flux of energy. Moreover, 

S l h(r,t) = (2£ + 3)& T0 y4 £] (r,f) + 

6ro{^[l [2l 4 _2l (r,f)} + b Tl — {£*[l [2] /jf 2 _2] (r,t)]} + 

TXl 

a L0 rnV ■ 4 +1] (r,f) + 

Oil— ^<r(M)[V/|£" 1] (r,t)] + 

7 / Z 

3(Z + 2)a T i-—/j^r? + 

TXl 

flfW), (33) 
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is a contribution which couples the thermal motion to the material motion (thermal-striction 
effect), R^\r,t) contains contribution of all the other higher-order fluxes £ + 2). The 
several kinetic coefficients are presented in Appendix B in Ref. [13]. 

Taking into account that 

|6rl1 = 5(l + x) |aTo1, (34) 

it follows that Maxwell times have the property that 

0h e+1 _ 5(T + 2)(1 + x) + £ 2 + 5 (£ + 1) , x 

6 ht ~ 5(£ + 3)(l + a;) + (^ + l) 2 + 5(£ + 2) ’ 1 ’ 

for £ = 0,1, 2, 3,..and where x = m/M. The ordered sequence 

Oho > Ohi > @h 2 > 0h 3 > ■ ■ ■ > @h( > @h e+ 1 > • • • (36) 


is verified in an analogous condition on the Maxwell times associated to the particle motion. 
It may be noticed that for large £, 9( goes as £ 2 , i.e., going to zero as £ —> oo. 

For the Brownian particle (x S> 1) the ratio in Eq. (35) tends asymptotically to (£+2 )/(£+ 
3) while for Lorentz particles (x <C 1) the ratio tends to [t' 2 + 5(2£ + 3)]/[(£ + l) 2 + 5(2£ + 5)]]. 
Considering that the Maxwell time of the first flux of mater, 9 ni is loqol^ 1 (cf. Eq. (22)) we 
can establish the relation between it and the Maxwell time for the energy 


0h o 


1 + x 
2x + 3 


'n\ i 


(37) 


what tell us that 0h o < 0 ni . For a Brownian particle 0h o = 0 ni /2 and for a Lorentz particle 
0h o = 0 n i/3. Moreover, for general £ we have that 


Oh, 


_ 5(1 + a;) _ „ 

5(1 + 2)(1 + x) + P + 5(1 + 1) 1 ’ 


(38) 


Therefore, recalling that the first flux of matter multiplied by the particle mass is the 
linear momentum density, as a general rule we can state that the energy relaxation time 
is smaller than momentum relaxation time. As already notice at the end of the previous 
Section, the hierarchy shown in Eq. (36) is quite relevant for the choose of the contracted 
description to be used. 

We proceed next to consider the case of a MHT at a quantum mechanical level for the 
cases of phonons and electrons in semiconductors. 
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IV. MAXWELL TIMES IN MHT OF PHONONS IN SEMICONDUCTORS 


We consider a system of longitudinal acoustic la phonons in a semiconductor in anhar- 
monic interaction with the accompanying transverse acoustic ta phonons. The sample is in 
contact with a thermostat at temperature T 0 . An external pumping source drives the LA 
phonon system out of equilibrium. The system is characterized at the microscopic level by 
the Hamiltonian 

H = Hos + Hob + Hsb + Hsp , (39) 

which consists of the Hamiltonian of the free LA phonons 

Hos = ^q(aqOq + 1/2), (40) 

q 

where is the frequency dispersion relation and the sum on q runs over the Brillouin zone, 
and the one of the TA phonons, which we call the thermal bath in which the LA phonons 
are embedded, given by 

H 0B = J2 m ^ b l b « + - 1 / 2 )’ ( 41 ) 

q 

where H q is the frequency dispersion relation; o q (o q ) and & q (6 q ) are the corresponding 
creation (annihilation) operators in mode q. 

Moreover, the interaction of the LA phonons with the thermal bath is given by 

Hsb = ^k.qOq/^+q^— k + he , (42) 

k.q 

where we have retained the only process that contributes to the kinetic equations 
(la^tta+ta) and we have neglected nonlinear contributions; M kq accounts for the cou¬ 
pling strengths. Finally Hsp is the interaction energy operator for the phonons and an 
external pumping source, to be specified in each case. 

At the macroscopic level (the nonequilibrium thermodynamic state) NESEF requires 
first to specify the basic variables that are to be used to characterize the nonequilibrium 
ensemble [7-10]. A priori, when the system is initially driven away from equilibrium, it is 
necessary to include all the observables of the system what is attained introducing many- 
particle dynamical operators [14,15]; in the present case it suffices to take only the single¬ 
phonon dynamical operators h q ,Q = ® q+ Q/ 2 a q-Q/2 in the second-quantization representation 
in reciprocal space. The two-phonon dynamical operator and higher-order ones can be 
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ignored because of Bogoliubov’s principle of correlation weakening [15]. Moreover, it would 
be necessary to also include the amplitudes cd and a q because their eigenstates are the 
coherent states [16], and the phonon pairs [17], both are disregarded because they are of no 
practical relevance for the problem considered here. Another basic microdynamical variable 
that needs be incorporated is the energy of the bath, and then the basic set is composed of 

{£q, £q,Q, H ob } , (43) 

with Q d 0, where we have separated out the so-called populations, z> q = a q a q , from those 
with Q d 0 which are related to the change in space of the populations (they are also called 
coherences [18]). 

Therefore, for the present case the nonequilibrium statistical operator is 

t 

g e (t) = exp{-S(t, 0) + J dt'e^'^-^Sif ,t'- t)}, (44) 

— OO 

where 

S(t, 0) = — ln£>(t, 0) = 

= 0(f) + F q(f) a K + F nQ( t ) a g+Q/2 a 9-Q/2 + PoHob , (45) 

q qQ^O 

S(t',t' — t) — exp — (t 1 — t)H^ S(t', 0) exp | — (t 1 — t)H^ , (46) 

Hqb is given in Eq. (41), 0(f) = lnZ(f) plays the role of a logarithm of a nonequilibrium 
partition function ^(f), and 0o = 1 /(^b7o). 

The average values of the microdynamical variables in set (43) over the nonequilibrium 
ensemble provide the variables that characterized the nonequilibrium macroscopic state of 
the system, which we indicate by 

Wf), ^qQ(f), E b}, (47) 

and the nonequilibrium thermodynamic variables conjugated to them are 

(F q (f), F qQ (f), A,}. (48) 

Going over to direct space (anti-transforming Fourier in variable Q) we obtain the space 
and crystalline momentum-dependent distribution function z/ q (r,f), in terms of which is 
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built the phonon Higher-Order Generalized Hydro-Thermodynamics which consists of two 
families of hydrodynamical variables, namely, the one associated to the quasi-particles (the 
phonons) motion (we call it the n-family) 

Mr, t), I„(r, t), {l}?( r, t)}} , (49) 

where £ = 2,3,..., and n(r,t) stands for the number of phonons at time t in position r, 
namely 

n ( r ,t) = ^2 iy n(r, t ), (50) 

q 

the hrst flux of this quantity 

In(r,t) = V q cu q z/ q (r,t) , (51) 

q 

and the higher-order fluxes 

4?( r >*) = X)“ M ( q ) I/ q( r ’ i ) ’ ( 52 ) 

q 

with £ = 2, 3,..., defining the Gorder flux (Grank tensor) where 

rt^(q) = [V q a; q : .. .£ — times ... : V q ca q ] , (53) 

is a Grank tensor consisting of the tensorial inner product of Gtirnes the group velocity of 
the q-mode phonon, VqCUq. 

On the other hand we do have the family associated to the energy motion (heat transport, 
and we call it the /?,-family) 

(h(r,f), I h (r ,t), (4 £] (r,f)}}, (54) 

where 

h( r, t) = ^ ftw q z/ q (r, t ), (55) 

q 

Ifc(r, t) = hWqVqWq i/ q (r, t), (56) 

q 

^(M) = q)t/ q (r,f) , (57) 

q 

which are, respectively, the energy density, its hrst (vectorial) flux, and the higher-order 
(£ — 2, 3,...) tensorial fluxes at time t in position r. 
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Consequently, the hydrodynamic equations of motion (evolution equations for the quan¬ 
tities above) are 

= J]/4 £1 (q)^z/ q (r,t) , (58) 

q 

where, for p = n, Kn\ q) = rt^(q) and for p = h, q) = q), recalling that £ = 0 

corresponds to the densities, £ = 1 to the vectorial fluxes, and £ ^ 2 for the higher-order 
fluxes. 

Evidently, all the evolution equations in Eq. (58) are determined by the unique equa¬ 
tion of motion for the single-phonon distribution function z/ q (r, t). It follows from the evo¬ 
lution equation which results in taking the average over the nonequilibrium ensemble of 
the quantum mechanical Heisenberg equation of motion for the microdynamical variable 
h q Q = a q+ Q/ 9 a q _Q/ 2 , (for practical convenience is calculated in reciprocal space), that is 

^q.qW = Tr j^[z> qQ , #]&(£) x g B | . (59) 

where p e {t) is the nonequilibrium statistical operator of Eq. (44), and qb is the canonical sta¬ 
tistical distribution of the thermal bath (the TA phonons) at the temperature T 0 . Performing 
such average is extremely difficult and then it is necessary to resort to the introduction of 
a practical nonlinear quantum kinetic theory [7-10,19,20], which, for the present case, can 
be used in the approximation consisting in retaining in the collision integral contributions 
up to second order in the interaction strength (Markovian approximation [9,19,21]). In that 
approximation it follows, after going over to direct space, a generalization of the so-called 
Peierls-Boltzmann equation given by [22] 

^z/q(r,f) = -V-Vq£h q z/q(r,f)-r q [z/ q (r,f)-z/°] 

+Jqs(r,t), (60) 


where J q s contains the effect of the presence of external sources, to be specified in each 
particular case, and 

u° = --- 

q g/3 0 fitJq _ 1 ’ 

is the distribution in equilibrium, 


4“ -fq ? 


( 61 ) 


with 




l + vi A + u, 


k+q 


f^k+q + fh 


Wn 


( 62 ) 
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which is the so-called self-energy correction, and 


r q — ^2 ^2 l^kqh(l + ^k‘ 4 + ,y k+q)^(^ k +q + ~ ^q) + T q , (63) 

k 

plays the role of the reciprocal of a relaxation time per mode, say r q 1 , (towards the equi¬ 
librium distribution), where 


'P-gasnv < 64 > 

is the distribution in equilibrium at temperature T 0 of the TA phonons. In this Eq. (63) we 
have written explicitly the known contribution due to the anharmonic interaction, and in a 
Mathicssen rule [23] are also indicated in T q the sum of the inverse of the relaxation times 
due to the presence of impurities, imperfections, boundary conditions, etc. 

Using Eq. (60) in Eq. (58), we obtain for the set of evolution equations that 


= 5^-^p l (q)[- V - V qWq]^q( r ,i) 
q 

-Eq Mr, t) - u°] + jjf(r, t ), (65) 


for both families of sets (49) and (54), p = n and p = h respectively. A closure for the 
set of these equations must be introduced, that is, to express z/ q (r, t) that appears on the 
right hand side in terms of the hydrodynamic variables in the sets of Eqs. (49) and (54). 
First it can be noticed that Eq. (65) involves an enormous set of coupled nonlinear integro- 
differential equations, he., both densities together with their fluxes of all orders are coupled 
through all these equations. To proceed further it is necessary to introduce a contraction 
of description , that is to say, to reduce the number of variables to be used, as described 
elsewhere [9,10,24], Moreover, the two families in the sets of Eqs. (49) and (54) are coupled 
by cross-terms that account for thermo-striction effects. In cases where these effects are 
not particularly relevant they can be disregarded and we obtain two independent sets of 
evolution equations, one for the n-family and the other for the /?,-family. 

Admitting that the just mentioned conditions are well satisfied, we introduce a contracted 
description of first order, and proceed to analyze the evolution equations for the energy 
density and its first flux. The resulting equations are [22] 

^-h(r, t) + V • I h (r, t) = a L (t)V • I h (r, t) 
at 

-d^(t)h(r, t ) + J hS ( r, t) , (66) 
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( 67 ) 


^I»(r,t) + V-^(r,t)=6PW-V/>(r,t) 

-6»f 1 (t)I A (r,t) + J IS {r,t ), 


where 



a L (t) = ^2 ^a q a/(q, t) ■ V q P q , 

q 

(68) 


bfit) = 22 (q, t) [V q Pq ) V q W q ] , 

q 

(69) 


a h (q,t) = ha; q z/q(t)[l + z/q(t)][^(hu;q) 2 z/q(t)[l + z/q(t)] , 

q 

(70) 


a/(q, t) = hw q Z/q(t)[l + Z/q(t)] ^^(hca q ) 2 Z/q(t)[l + Z/q(t)] X 

q 



1 _1 

[VqCJq * Vq^q] VqCJq , 

(71) 

P q and r q are given in Eqs. (62) and (63), and 



e F q (t) _ 1 ’ 

(72) 

with P q of Eq. 

(48), we recall that 



4 2] (Q,t) = X] ^q[V q W q : VqWq]z/q Q (t) , 

(73) 

and Jhs and 3 is are the contributions due to the coupling with the external sources 

that are 

present which 

are to be specihed in each case being considered. We recall that [• ■ 

. : ...] 

stands for inner tensorial product. 


Finally, 

9 h 1 ( t ) = 5^">h( q ,*) r q* 

q 

(74) 

and 

Oj 1 ^) = ^W/(q,f) r q, 

a 

(75) 

with 

n 

^(q.t) = v ^ ) > (1 n + ; J r 

2^(PUq/)V q /(l + Z/q/) 

q' 

(76) 


... 7 n ,1 (^q)V q (l + Pj |VqWq| 2 

/lq ’ ^ E(^q') 2 ^(l+^q')|V q ^| 2,: 

(77) 


q' 
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which have the property of normalization, that is, it is verified that 


J2w h (q,t) = l, (78) 

q 

and 

5^wj(q,*) = 1 , (79) 

q 

Equations (74) and (75) define the reciprocal of Maxwell times associated to energy and 
its first flux, respectively, in this MHT of order 1. We can see that each Maxwell time follows 
a kind of Mathiessen rule: its inverse is a superposition of the reciprocal of the relaxation 
times of each mode r" 1 = T q , multiplied by a weighting function w(q,t). Because of 
the normalization condition (Eqs. (78) and (79)) it can be considered each w(q, t) as the 
probability at each time t of the contribution of the mode-relaxation time r q . 

In a Debye model both Maxwell times are equal, i.e., 


o- h \t) = ej\t) 


^(sg) 2 Z/q(l + Z/ q )r q 

EN^qfi+b) 


(80) 


and we recall that q < qoebye- 

Taking into account that |V q o; q | ~ s for small values of q and that |V q o; q | < s for 
intermediate to large values of q, we can estimate that in general 6h > 0j. It is conjectured 
that we would do have a hierarchy 


> &i > > • • • ) 


quite similar to the case of the classical fluid of previous sections. 
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V. MAXWELL TIMES IN MHT OF ELECTRONS IN SEMICONDUCTORS 


We consider the case of a n-doped polar semiconductor, e.g., n-GaAs or n-GaN, whose 
Hamiltonian quantum mechanical operator is 

H = Hoe + H ep h + Ho p h + H ep , (81) 

where Hq c is the Hamiltonian of the carriers (a concentration n of electrons in conduc¬ 
tion band) taken in the effective mass approximation, with energy dispersion relation 
e k = h 2 k 2 /2m*. H e ph is the Hamiltonian of the interactions of carriers and optical and acous¬ 
tical phonons: we consider, in this case of polar semiconductors only the polar (Frohlich) 
interaction with LO phonons, which is by large the predominant process of relaxation of the 
carriers to the lattice. H 0p h is the Hamiltonian of the free LO phonons. H eip accounts for the 
interaction of the carriers with external pumping sources. We assume that the LO phonons 
remain in equilibrium with an external reservoir at temperature T 0 . The polarization effects 
(due the Coulomb interaction) modifying the mean held potential of the band states is not 
included, that is, we disregard the contribution from the plasmon states. 

Indicating by c k (cj <; ) the creation and annihilation operators in band states k (we omit 
the spin index), as in the case of phonons of the previous Section we introduce as basic 
microdynamical variables for the corresponding NESEF the quantities 

{'h k Q = c k +Q/2 C k-Q/2; Hoph} ■ (82) 

For <5 = 0, h k = <4.c k is the operator occupation number in state k, and the others with 
Q ^ 0 accounts for variations in space. For the LO phonons we have taken Ho p h once, as 
said, we consider them in constant equilibrium at temperature To and therefore described 
by a canonical distribution. Hence, the statistical operator in NESEF of the electrons is 

t. 

g e (t) = exp | — S(t, 0) + J dt'e 6 ^ S(t',t'— t) j, (83) 

— OO 

where 


S(t u t 2 ) = lng(t 1 ,t 2 ) = 

— 0(G) + y> k (tiK(t 2 ) + ^ F kQ (ti)n kQ (t 2 )], (84) 

k Q^O 
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where we recall that 0 ensures the normalization of both g and g e , and F k and FkQ are the 
nonequilibrium thermodynamic variables conjugated to the basic ones h k and h k Q (Q / 0). 
In Eq. (84) ti refers to the evolution in time of the nonequilibrium thermodynamic variables 
and £2 to the evolution in time of the dynamical variables (Heisenberg representation). 

Following the formalism we proceed to the derivation of the kinetic equation for the 
variable 

«kQ(£) = Tr{h kQ £ e (£). (85) 

In the Markovian approximation after a calculation in reciprocal space and next going 
over to direct space it follows for 

/k(r, t) = Y n kQ(t)e* Q ' r , (86) 

Q 

that 

d x 

T^/k(r, t) + ^V k e k • V r / k (r, t) = 

O-tt 

-jr Y l C q \ 2 {K(t) + l]/k+q(r,£)[l - / k (r, t)} - 

q 

^ q (t)/k(r,£)[l - /k+ q (r,£)]}h(e k+q - e k - £ku q ) - 

-f~Y l C q| 2 {KW + ii/kMMl - /k- q (r ,t)] - 

q 

^ q (t)/k-q(r,£)[l - / k (r,£)]}5(e k _ q - e k + hca q )} + 

^•k(r,£), (87) 

where C q is the matrix element of the Frohlich clectron-LO phonon interaction, ta q the 
frequency dispersion relation of the LO phonons, we recall that e k = fr 2 k 2 /2m*, Jsk the 
interaction with external sources to be specified in each case, and z/ q is the population of 
the LO phonons in equilibrium at temperature T 0 . 

As in previous Sections we introduce the hydro-thermodynamic variables 


i 1 p ] (r, t) = Y K p ] ( k )/k( r > *), 

k 

(88) 

k ) = ~ .. .£ — times , 

m* m* 

(89) 


(90) 


19 



corresponding to the families of particle motion (p = n) and energy motion (p = h). The 
evolution equations are evidently given by 

|^ ] ( r ,f) = J2 K p ] (k)^f k (r,t ), (91) 

k 

In what follows, for better visualization, we restrict the analysis to a description of order 
2 of the heat transport introducing the reduced set of basic variables consisting of 

{h{ r, t); 4(r, t); I h2 ( r, t)} , (92) 

where 

h 2 k 2 

M r 4) = E 2^k(r,t), ( 93 ) 

k 

__ fcb 

k(r,t) = / k (r,t), (94) 

and we separate the second order flux /j^(r,f) in the trace part and the remaining traceless 
part which is neglected. The trace of 4 2 1 is 

4 2 (M) = /k( r 4) • (95) 

k ' ' 

The auxiliary statistical operator is then given by 
g(t, 0) = exp | - <j>(t, 0) 

- 5^[-PX(Q, t)h( Q) + y? h (Q, t) • 4(Q, t) 

Q 

+4i2(Q4)4 2 (Q)]| , (96) 

which follows from Eq. (84) after taking, to be consistent with the contraction introduced, 
that 

FkqW = ^ Fhl ' Q ’ t) + ^ hk ' vM ’ t) + 

/ fjU \ 2 

— F h2 (Q,t). (97) 

\ m*) 

The equation of evolution for these hydro-thermodynamic variables in the absence of 
external sources are 

c) 

—h(r, t) = -V ■ 4(r, t) - 9^h( r, t) - a 0 4 2 (r, t ), (98) 
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(99) 


d_ 

dt 


Ift(M) 


-vi h 2 {v,t)-e^i h {T,t), 


—4 2 (r, t) = ^ ( r , t) - a 2 h( r, t) - 9j^I h2 ( r, t). (100) 

The expression for the trace of the divergence of the third order flux in terms of the basic 
macrovariables, namely 


Tr {V ■ 4 3, (r,f)} = JjV-In(r ,t) -a 2 h(i,t), 


( 101 ) 


was obtained resorting to Heims-Jaynes perturbation expansion for averages [25] around the 
homogeneous state where 


a 2 


1 

m*/36 l2 ’ 


( 102 ) 


after using that 


/k = 


exp {(3 h 2 k 2 /2m*} + 1 ’ 
1 




k R T* ’ 


(103) 

(104) 


with T* being the nonequilibrium temperature (quasitemperature) of the electrons, and in 
Eq. (84) making the identification 


F k = 


h 2 k 2 

2m* 


p(t), 


(105) 


In Eqs. (98) to (100) are present the Maxwell times 9 h , 9 Ih and 9 Ih2 , associated to the 
energy, its vectorial (first) flux (current of heat) and the second flux respectively. They are 
given by 


where 


and 


q-- 1 - —fir' 1 

T h. — 5 1 

(106) 

A, r -1 - — 

(107) 

—e 2 ^ 1 ~ — — f 3 ^ 1 

28^ h 25 5 ^ 0 ’ 

(108) 

A k B T* ’ 

(109) 

To 1 = i 

(110) 


in the calculation being used the parameters corresponding to GaAs given in Table I. 
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TABLE I: Parameters characteristic of GaAs 


Parameter 

Value 

electron effective mass, m* 

0.067too 

static dielectric constant, £o 

12.91 

optical dielectric constant, £oo 

10.91 

LO phonon energy, Hulo 

37.0 meV 

mass density, p 

5.31 g/crn 3 


Accordingly, it follows that the comparison of the three values follows from 


Th 




Th 

ZIh 

Tl 


875 k B T* 
14 huj 0 


k B T* 

hu 0 

56 k B T 


( 111 ) 


( 112 ) 


(113) 


125 Tiuiq ’ 

and, therefore, their hierarchy depends on the ratio hoJo/k B T *, that is, on the ratio of the 
energy of LO-phonon modes to the nonequilibrium thermal kinetic energy of the carriers. 
Numerical values of the three Maxwell times for several values of the carriers’ quasitemper¬ 
ature are given in Table II. 


TABLE II: Numerical Values of Maxwell Times for GaAs 


T* (K) e h (ps) 

e i h (Ps) 

e i h 2 (P S ) 

100 

23.7 

33.76 

193 

150 

3 

6.4 

36.5 

200 

1.24 

3.5 

20 

300 

0.57 

2.45 

19 

500 

0.43 

3.06 

17.5 

1000 

0.43 

6.12 

3.5 


It can be noticed that these Maxwell times are in order of picoseconds to femtoseconds. 
We stress that T* is the carriers’ quasitemperature, that is, the nonequilibrium temperature 
which is a measure in Kelvin degrees of their nonequilibrium energy [26]. 
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Inspection of Table II tells us that it is verified that Oh < 6j h < 0j h21 with this ordering 
being inverted for values of £ sufficiently larger than 1, i.e., for carriers’ quasitemperatures 
sufficiently larger than Einstein temperature (T E = hujo/ks) of the LO phonons, roughly 
T* > 432 Kelvin in GaAs. 
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VI. CONCLUDING REMARKS 


In the framework of Mesoscopic Hydro-Thermodynamics applied to the study of the 
hydrodynamic motion of a classical molecular fluid, and of the ”fluids” of phonons and 
carriers in intrinsic and doped semiconductors, we have characterized the so-called Maxwell 
times, or, better to say, generalizations of Maxwell original proposal [1] as described in the 
Introduction. These Maxwell times are associated with the fluxes of all orders present in 
MHT (fluxes of particles - or quasi-particles - and of heat). A complete characterization is 
obtained in the case of a classical fluid (of molecules, polymers, etc.), presented in Section II 
and III. It follows the quite important result that they satisfy a hierarchy showing that they 
are increasingly diminishing as the order of the fluxes increases. This allows for establishing 
criteria for the choice of the contraction of description to be used in each case. 

In other words to introduce an appropriate - for each case - contraction of description: 
this contraction implies in retaining the mformation considered as relevant for the problem 
in hands, and to disregard irrelevant information [27]. 

Elsewhere [24] it has been discussed the question of the contraction of description (re¬ 
duction of the dimensions of the nonequilibrium thermodynamic space of states), where a 
criterion for justifying the different levels of contraction is derived: It depends on the range of 
wavelengths and frequencies which are relevant for the characterization, in terms of normal 
modes, of the hydro-t her mo dynamic motion in the nonequilibrium open system. It can be 
shown that the truncation criterion rests on the characteristics of the hydrodynamic motion 
that develops under the given experimental procedure. 

Inclusion of higher and higher-order fluxes implies in describing a motion involving in¬ 
creasing Knudsen numbers per hydrodynamic mode, that is, governed by smaller and smaller 
wavelengths - larger and larger wavenumbers - accompanied by higher and higher frequen¬ 
cies. In a qualitative manner, we can say that, as a general “thumb rule”, the criterion 
indicates that a more and more restricted contraction can be used when larger and larger 
are the prevalent wavelengths in the motion. Therefore, in simpler words, when the motion 
becomes more and more smooth in space and time the more reduced can be the dimension 
of the space of basic macrovariables to be used for the description of the nonequilibrium 
thermodynamic state of the system. It can be conjectured a general criterion for contrac¬ 
tion, namely, a contraction of order r (meaning keeping the densities and their fluxes up to 
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order r) can be introduced once we can show that in the spectrum of wavelengths, which 
characterize the motion, predominate those larger than a “frontier” one, X lr,r+1) = V 2 9 r 0r +1 
where v is of the order of the thermal velocity and 9 r and 9 r+ i the corresponding Maxwell 
times associated to the r and r + 1 order fluxes. 

Section IV has been devoted to the characterization of Maxwell times in the thermal 
motion of phonons in intrinsic semiconductors. It is shown the interesting results that the 
Maxwell times have an expression for their inverse in the form of a weighted Mathiessen-like 
rule involving the inverse of the relaxation time of all phonon modes. We can say that it is 
expected once the hydrodynamic motion is a collective one composed of the contribution of 
all the phonon modes. As in the previous case Maxwell times follow a hierarchy involving 
ever decreasing values as the order of the fluxes increases. 

Finally, in Section V, for the case of doped polar semiconductors, the Maxwell times 
associated to the hydrodynamic motion of the electrons in Bloch bands have been obtained. 
They are mainly determined by the presence of the polar Frohlich interaction between elec¬ 
trons and phonons and they follow a hierarchy of values, which is dependent on the ratio of 
the phonon energy Hujq with the electron thermal (kinetic) energy: it implies in a series of 
decreasing values as long as the quasitemperature of the “hot” electrons is higher than the 
Einstein temperature of the phonons. 
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